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The first observations of the rare decays Ks — ► vr e + e and Jfs — > 7r /^ + /^ have been made 
i-^ by the NA48/1 collaboration at the CERN SPS accelerator. From high intensity Ks data 

collected during the 2002 run, clean signals of 7 Ks — » 7r°e + e _ events and 6 isTs — » 7r /j + ^ _ 
events were observed, giving branching ratio measurements of BR(Ks — > 7r°e + e~) = 5.8JI24 x 
10~ 9 and BR(K S -> n°^ + u') = 2.9±J;| x 10" 9 . These results constrain the indirect CP 
violating component of the corresponding i^t decays. Other recent results from NA48 are 
also presented. 



1 Introduction 

The NA48 experiment was originally designed to measure the CP violation parameter Re(— ) 
via a high statistics comparison of Ks — ► vt7t and Kl — ► vr-7r decays. A beamline consisting 
of simultaneous Iv^ and iv^, beams was used, with a tagging system allowing discrimination 
between the two types of kaon. Rare Kt, decay and precision measurements could be carried 
out in parallel. In 2002, for the NA48/1 phase of the experiment, the K l beam was removed 
and the proton intensity on the Ks target was increased by a factor ~ 1000, allowing very rare 
Ks and hyperon decay searches to be pursued. 

The kaon beams were produced using 400 GeV/c protons provided by the Super Proton 
Synchrotron (SPS) at CERN, incident on two separate beryllium targets. The NA48 detec- 
tor consisted of the following principal sub-detectors: a magnetic spectrometer consisting of 4 
drift chambers separated by a dipole magnet, a high resolution liquid krypton electromagnetic 
calorimeter, an iron scintillator hadronic calorimeter and a muon system consisting of 3 planes 
of scintillator shielded by 80cm thick iron walls. 1 
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Figure 1: (a) m ee distribution of the Dalitz decay and conversion background. Reasonable agreement between 
data (points) and Monte Carlo (curves) can be seen, (b) 2001 data showing the distribution of the Kl — > ee-yy 

background across the signal region. 

2 Physics Motivation for Ks — ► tt°1 + 1~ Searches 

The decay mode Kl — » ts^vv is direct CP violating and can be used to determine the parameter 
7] in the Wolfenstein parameterisation of the CKM matrix. This channel is theoretically very 
clean but is experimentally challenging due to the missing energy from the neutrinos. A second 
possibility to measure rj is from the decay Kl — > ir°l + l~. This decay is easier to measure 
experimentally as all the decay products are detectable, but is more complicated theoretically as 
it has contributions from CP conserving and both direct and indirect CP violating components, 
which also interfere. The CP conserving component can be predicted from a measurement of 
Kl — »• 7r°77 while the indirect CP violating component can be predicted from a measurement 
of Ks — > 7r°Z + Z - . Any measurement of BR(Kl —> tt°1 + 1~), and consequently, 77, requires the 
measurement of BR(K$ — ► 7r°Z + Z — ) to disentangle the direct from the indirect CP violating 
components. 



3 Analysis Strategy for Ks 



7v°e+e- 



and Ks 



7rVV~ 



Both Ks — ► 7T e + e~ and Ks — » 71 fx + \x~ were predicted to have very small branching ratios, 
(~ 1CP 9 ) and so, even with a flux of ~ 3 x 10 10 Ks decays, only a handful of events were 
expected from either channel. Therefore, to avoid biasing the results, a blind analysis procedure 
was employed. A 2.5a mK x 2.5a n o signal region and a 6.0a mK x 6.0a n o control region were 
defined. Both regions were kept masked during the background studies. Only after the back- 
ground contributions to the signal had been estimated and the cuts fixed was the control region 
unmasked. Final selection changes could then be made if required before the signal region was 
unmasked. 



4 K s -*■ 7r°e+e- 

4-1 Signal Selection 



Ks — ► 7r u e + e~ candidates were selected with 40 < E& < 240GeV and within 2.5 Ks lifetimes 
of the end of the final collimator. Two oppositely charged tracks, with E/p > 0.95 and forming 
a good vertex, were required, as well as two clusters in the calorimeter that were not associated 
to tracks. 



0) 

O 



0.6 






0.55 












0.5 


gjj#- 






" 




0.45 






0.4 


1 , , , , 1 , 





0.1 



0.15 0.2 



0.25 



0.3 0.35 

m ee /GeV 



> 

CD 

o 



0.65 
0.6 

0.55 
0.5 

0.45 
0.4 

0.35 



- 


0.5 






- 


0.49 






0.134 


0.136 


- 




















0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 

m/GeV 

yy 



Figure 2: The (a) m 11 vs. rn ee and (b) m 11 vs. uik mass planes showing Ks — > 7r°e + e candidate events. The 
large background from dalitz decays and photon conversions can be seen at low m ee 



.2 Background 



Backgrounds for the Kg — > tt e + e~ decay were divided into two categories: physical background 
(arising from a single kaon decay), and accidental background (arising from two separate, over- 
lapping kaon decays). Three significant sources of physical background were identified. The first 
came from Dalitz decays (ir° — ► ee^/) and photon conversion in the various parts of the detector. 
These backgrounds were studied extensively using Monte Carlo simulation (see fig. la). Reason- 
able agreement with data was found in the m ee distribution around the 7r° mass. A conservative 
cut of m ee > 0.165 GeV/c 2 was applied to remove this background, with a corresponding 48% 
loss in acceptance. 

The second significant physical background was identified as Kl — > ee'yy. This background 
spread across the m 77 signal region as can be seen in fig. lb. Using Kl data taken in 2001, the 
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background was estimated to be 0.08_ 002 events. 

The final physical background came from various hyperon decay channels, 
contained a significant (~ 10 9 ) flux of neutral hyperons and the channels S° — 
H° — > A(pe~ v)tt° and H° — > S + (p7r°)e _ ^ were identified as possible background sources. A cut 
on the momentum asymmetry of the decay products reduced these backgrounds to a negligible 
level. 

The accidental background was dominated by overlapping fragments of two decays, e.g 
(Kl — ► -K^e^v) + (Ks — ► 7r°7r°). This background was suppressed by cutting on the spread 
of track and cluster times (|Ai| < 3ns). To estimate the remaining accidental background, the 
timing cuts were relaxed and the time sidebands investigated. No events were found in the 
out of time (3ns < |Ai| < 50ns) signal region. Therefore, an extrapolation was made from the 
out of time control region to the out of time signal region, taking into account the shape of the 
background from Monte Carlo simulation. A further time extrapolation could then be performed 
to get the final background estimate for the in time region (|At| < 3ns) of 0.069 events. 



4.3 Result 
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The overall background estimate was 0.15_ ' 04 events. The control region was unmasked and 
no events were found. No changes to the selection were needed before the signal region was 
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Figure 3: (a) The kaon mass region with the Kl — > ty°tv + tv background. Reasonable agreement between data 
and monte carlo (24 times the data) is shown, (b) Distribution of out of time events in the ttik vs. m 77 plane. 



unmasked. A signal of 7 events was found, giving the result: 2 

BR{K S -» vr°e + e-) = (5.8±|;|(stot) ± 0.8(syst)) x 1(T 9 



(1) 



where the systematic error is dominated by uncertainties arising from the extrapolation to the 
full m ee region. The candidate events are shown in fig. 2. 



5 K s — ► 7rV+/x- 

5.1 Signal Selection 



Kg — > 7i n + n~ candidates were selected with 60 < Ej( < 200GeV and within 3 Kg lifetimes of 
the end of the final collimator. Two oppositely charged tracks were required that formed a good 
vertex, as well as two clusters in the calorimeter that were not associated to tracks. The tracks 
had to leave minimum energy in the calorimetry and have associated hits in the muon detector. 



5. 2 Background 

As with Kg — > 7r°e + e _ , the backgrounds were divided into physical and accidental sources. 
Three dominant sources of physical background were identified. The first was Kl — ► 7r°7r + 7r~ 
where the two charged pions had decayed in flight. This background was studied extensively in 
Monte Carlo, (see fig. 3a) and estimated to be < 0.019 events. The second significant background 
was Kl — ► /U + /i~77. This, again, was studied with Monte Carlo and estimated to contribute 
0.04 ± 0.04 events to the overall background. The last physical background came from hyperon 
decays. As with Ks — > 7r°e + e _ , a cut on the momentum asymmetry of the decay products was 
used to reduce this background to a negligible level. 

The accidental background was treated similarly to K$ — > 7r°e + e _ , by loosening the timing 
cuts and extrapolating into the in time (|At| < 1.5ns) region. 6 events were found in the out of 
time (—115ns < At < 60ns) signal region, leading to a background estimate of 0.18_ ' 11 events. 
The distribution of these events in the mass plane is shown in fig. 3b. 



5.3 Result 
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The overall background estimate was 0.22_ ' 12 events. The control region was unmasked and 
no events were found. No changes to the selection were needed, the signal region was unmasked 
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Figure 4: The (a) m 77 vs. m ee and (b) m 77 vs. ttik mass planes showing Ks — > 7r°/i + /x _ candidate events. The 
events at large m 77 are consistent with ift, — ► 7r°7r + 7r - background and the single event at low m 77 is consistent 

with accidental background 



and 6 events were found. This led to the result: 

BR(K S -> 7rV + M~) = (2.9iJ:|(fltat) ± 0.2(syst)) x 10~ 9 (2) 

where the systematic error comes from uncertainties in the normalisation, obtained from Ks — ► 
7r + 7r _ decays. The candidate events are shown in fig. 4. 



6 Interpretation of the Ks — ► 7r°l + l Measurements 

The form factor for Ks — ► ir°l + l~ decays can be described in Chiral Perturbation Theory by a 
first order polynomial: 

W(z)~G F m 2 K (a s + b s z) (3) 

where z = mf l /m? K . This leads to the following predictions for the Ks — ► tt°1 + 1~ branching 
ratios? 

BR(K S -> vr°e + e") = [0.01 - 0.76a 5 - 0.216 5 + 46.5a| + 12.9a s b s + 1.44&|] x 10~ 10 (4) 

BR(K S -> ttV + /"~) = [0.07 - 4.52a 5 - 1.506s + 98.7a| + 57.7a 5 6 5 + 8.956|] x 10~ n (5) 

The Vector Meson Dominance (VMD) model predicts the ratio bs/as = 0.4, allowing \as\ to be 
extracted for both channels: 



BR(K S -» vr°e + e _ ) ~ 5.2 x 10~ 9 a| 
BR(K S -» ttVV - ) - 1.2 x 10~ 9 a| 



1.06+^? ± 0.07 



|asUo^ = 1.55i°;l±0.05 



(6) 
(7) 



These results for |as| agree within errors. 

By combining the ee and n^i results in a log- likelihood fit, as and bs can be determined 
separately. As can be seen from fig. 5, the observed Ks — ► 7r e + e~ and Kg — > 7r°// + /i _ rates are 
compatible with both each other and the VMD model. 





Figure 5: (a) Allowed regions of as and bs determined from the observed number of Ks —* 7r°e + e~ and Ks — > 
Tr°/j, + fjT events separately. The region between the inner and outer elliptical contours is the allowed region at 
68% CL. (b) Allowed regions of as and bs for the Ks — > 7r°e + e~ and Ks — ► 7r () /i + ^~ channels combined. The 
contours delimit the la and 2a allowed regions from the combined log-likelihood. The dashed straight line in 
both plots corresponds to bs = OAas, as predicted by the VMD model. 



7 Implications for Kl — ► 7r°Z + Z 

The CP violating component of the Kl — > tt°1 + 1~ branching ratio is given by: 

BR(K L -» TT°l + r) C pv x 10 12 
where Com is the direct CPV component, Cind 



Cind ± Cmix ( in -4 / 






(8) 



£-R(tf s -> 7r°/+r) is the indirect CPV 
component, Cmix ~ \JBR(K$ — > 7r u /+/ _ ) is the interference term and Im(At) = Im(V^*Vfd). 

The indirect CPV component and the interference term depend on the parameter \as\- Using 
the measured values of \as\ given above and taking 4 Im(At) = (1.36 ± 0.12) x 10 -4 gives the 
following predictions for the central value of the Kl — > tt 1 + 1~ branching ratios: 
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BR(K L ->1C e )cPV x 10 ±Z ~ 17i n dircct ± 9 int erference + 5<lirect 
BR(K S — > VT H + fJ>~)cPV X 10 ~ 9 ind i rcct ± 3i ntcr f C rcncc + Idirect 

where the sign ambiguity reflects the uncertainty in the sign of as- 



(9) 
(10) 
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it e^wy) comes from the KTeV collabora- 
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v) = (0.908 ± 0.008 



+0.013^ 
0.012 J 



(11) 



The most precise published measurement of BR(K° 
tion? 

BR(K° -» TT ± e T u- f )/BR(K° 

which is significantly below current theroretical predictions. This branching ratio has now been 
measured by NA48 using data from 1999. In order to reduce the uncertainties in the measure- 
ment due to radiative corrections, a model independent method was developed. The PHOTOS 
6 package was used for the initial Monte Carlo generation. The Monte Carlo events were then 
weighted with the 9* distribution, (the angle between the electron and the photon in the kaon 
rest frame). This gave good agreement between data and Monte Carlo for all distributions. The 
preliminary branching ratio measurement using this model independent method was found to 
be: 



BR(K° -» w ± e^uj)/BR(K° 



7 r ± e q F 



v) = (0.964 ± 0.008 



+0.012\ 
-O.OllJ 



(12) 



NA48 
(MI MC) 

KTEV 
NA31 



FFS 
Donsel 
ChPT 0(p 4 ) 



0.875 0.9 0.925 0.95 0.975 1 1.025 1.05 1.075 
BR(K e3y , 20°, 30 MeV)/BR(K e3 ) 



Figure 6: The measured value of BR(K° — > -k e^ vy) / B R(K {) — » ty e^v) compared with the published result 

from the KTeV collaboration and theoretical predictions. 



This is in good agreement with recent theoretical predictions, (Fearing, Fischbach and Smith 
(FFS) 7 ^, DonceP and Chiral Perturbation Theory calculations 10 ' 11 ), as can be seen from fig. 6. 
If the FFS method is used to apply the radiative corrections, as in the KTeV analysis, the result 
is found to be in agreement with the KTeV measurement. 
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A high precision measurement of the branching ratio and form factors for the Kl 
decay was made using Kl data from 2001. This decay provides a good test of Chiral Pertur- 
bation Theory (CHPT) predictions for long distance meson interactions and the form factor 
measurements allow the determination of the CHPT parameter, L3 . 12,13 

The main background to this decay was Kl — ► ir°ir + Tr~ where a pion was misidentified as 
an electron. To minimise this background, a neural network was used to distinguish between 
pions and electrons that took into account geometric characteristics of the showers and tracks. 
This reduced the background to 1.13%. 

A signal of 5464 events was found with a background of 62 events. This gave a preliminary 
branching ratio measurement of: 



BR(K L -> Tr^e*!/) = (5.21 ± 0.07 stat ± 0.09 SJ/st ) x 10" 



(13) 



Assuming a V-A structure for the matrix element, the form factor parameters were then mea- 
sured to be: 



f s = 0.052 ± 0.006 stat ± 0.002 



syst 



fp 



-0.051 i0.011.tat ±0.005 



syst 



X q = 0.087 i 0.019 stat i 0.006 



syst 



h 



-0.32 i0.12 sta ti 0.07, 



1 syst 



This led to the following value for L3: 

L 3 = (-4.1 ±0.2) x 1(T 3 (14) 

10 S — ► A7 

Using data from a short test run from 1999 with a high intensity Ks beam alone, 730 H — > 
A7 events were found with a background of 58.2 ± 7.8 events. This led to a branching ratio 
measurement of: 14 

BR(~ -> A 7 ) = (1.16 ± 0.05 stat ± 0M syst ) x 10" 3 (15) 

The statistics were large enough to enable a significant measurement of the decay asymmetry 
to be made: 

a(H -► A 7 ) = -0.78 ± 0.18 stai ± 0.06 SJ/st (16) 

This is the first evidence for a non-zero asymmetry in this decay mode. A much larger sample 
from the 2002 run is currently under analysis. 
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